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I present recent results related to the measurement of the top quark mass, using pp collisions recorded with the 
D0 Detector at the Fcrmilab Tevatron Collider. The results are: A direct measurement of the mass difference 
between top and anti-top quarks, Measurement of the top quark mass in the lepton+jets channel and in the 
dilepton channel. 



1. Introduction 



The D0 detector, at the Tevatron, is a gen- 
eral purpose particle detector, composed of tracking, 
calorimeter, and muon systems. The Tevatron collides 
opposing beams of protons and anti-protons, each 
with an energy of 0.98 TcV, the 1.96 TeV center of 
mass energy an energy makes the Tevatron the first, 
and still only, top quark factory. 

Top quarks at the tevatron are produced both singly 
and in pairs. However the single top production mode 
has both a smaller cross section, and more importantly 
a large irreducible background. When produced in 
pairs, both top and anti-top decay into a W boson and 
b quark, the W boson can then decay into a lepton and 
neutrino or a pair of quarks, and finally the b quarks, 
and any light quarks (u,d,s) produced by the W decay 
undergo particle showering and hadronization. 

Since 2001 the Tevatron has produced - 40000 tt 
pairs, though only a fraction of those events are useful 
for the mass analyses. 



2. Lepton+jets Event Selection 



There are two lepton+jets channels used for the 
mass analysis, electron+jets (e+jets) and muon+jets 
(/i+jets). Events in both channels are required to 
have a single isolated lepton with pt > 20 GeV and 
\rj\ < 1.1 (|fy| < 2 for muons), to have four jets with 
Pt > 20 GeV and \t)\ < 2.5, and to have large missing 
transverse energy JS T > 20 GeV (fi T > 25 GeV for 
muons). 

The main backgrounds in the lepton+jets channels 
come from multijet QCD, and VT+jets production. 
The QCD signal is reduced by placing quality cuts on 
the leptons, and also by placing a cut on the azimuthal 
angle between the lepton and the fix- The H^+jets 
background is reduced by requiring that one jet be 
tagged using a neural network discriminant trained to 
distinguish b-jets from light-jets. 



3. Top Mass Difference Measurement 

At the Tevatron, we can measure the top mass to a 
high precision, using three different channels, all-jets, 
lepton+jets, and dileptons. In all cases we assume 
that the mass of the top and anti-top are the same, 
such an equality is demanded by CPT (charge, parity, 
time) invariance [![. However, we know that studies of 
C and P individually led to the discovery that neither 
symmetry is conserved by the weak interaction, and 
similarly a violation of CPT would require new physics 
to explain. 



3.1. Matrix Element Method 

The top mass measurement uses the matrix element 
method (ME). The ME uses probability densities de- 
termined by integrating over the cross section for both 
the signal (tt), and the main background (W+jets). 
A likelihood is constructed having the form: 



£(|u,$ri! {u) 



(1) 



N evt 



JJ [ftPsig(x l ;mt,m i ) + (1 - ft)Pbk a (%i)] 

i=l 

Where xi represents the lepton and jet momenta 
measured in the detector and ft is the fraction of tt 
events in the sample. The signal probability P S i g has 
the form: 



P S ig(x;m t ,mi) = Acc(x) x 



(2) 



x / d n o(y;mt,mi)dq 1 dq 2 x f(qi)f(q2)W(x,y) 
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(2tt) 4 |.M(H£E -> Ud 



(3) 

t;$u4d)l e 



2qiq 2 s 



Where Acc(x) corrects for the acceptance of the 
event selection, a is the total cross section, y repre- 
sents the momenta of the tt decay products, \M\ e is 
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the matrix element of the process qq — > tt, d$>e is the 
6-body phase space integral, q\ and q2 are the mo- 
menta of the incoming partons, and W{x, y) is the 
transfer function which is the probability of parton 
level momenta y producing the measured momenta x. 
The background probability (Pbkg) has a similar form. 

For convenience we transform the original vari- 
ables (mt,mf) are transformed into two new variables: 
Am t = (m t — mj) and m sum = (m t + mj)/2. When 
quoting a final value for Am t , m sum is integrated out, 
with the final likelihood having the form: 



remove any bias on the measurement. The pseudo- 
experiments are formed using simulated events, and 
contain signal and background events consistent with 
that expected in data. The measurement procedure 
is performed for each pseudoexperiment in the same 
way as for data. 

The average of Am™ easmcd is compared to Am* rue 
and the small bias on A?7i™ easured [ s calibrated out 
(see Fig. [2]). Similarly, we look at the RMS of the pull 
{PULL = (Amf easured - Am* ruo )/<5(Am t ) mcasured ), 
and the bias on the measured error on A m m casured is 
removed (see Fig. [2]). 



£«M = J rWC(-$u,W ( 4 ) 

3.2. Simulation of Signal and 
Background Samples 

Signal samples for this analysis have been generated 
using the PYTHIA Monte Carlo event generator 
Top and anti-top masses are given values of 165, 170, 
175, and 180 GeV, with 14 samples in total (see Fig. 
P. 

The W+jets background was generated using the 
leading order matrix element generator ALPGEN 3], 
with PYTHIA used to shower the partons generated 
by ALPGEN. 
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Figure 1: mt and mj values for simulated signal. 
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Figure 2: < Amr asurcd > and pull RMS vs Am? 



3.4. Systematic Uncertainties 

Systematic uncertainties are determined by vary- 
ing the composition of the pseudoexperiments. The 
largest systematic errors come from signal modelling, 
QCD, and jet resolution. The total systematic error 
of 1.15 GeV is half the size of the measured statistical 
error 3.44 GeV, which is to say that this measurement 
is statistics limited. Table Q] contains a list of system- 
atic uncertainties. 



3.5. Measurement in Data 



The measurement in data is 



e + jets Am t = 0.33 ± 5.03(stat) GeV (5) 
/i + jets Am t = 6.74 ± 4.71 (stat) GeV 

comb Am t = 3.75 ± 3.U(stat) ± I.15(syst) GeV 



3.3. Performance on Simulated Data, 
Calibration of Method 

We have constructed ensembles of pseudoexperi- 
ments to test the performance of the method, and 



Figure [3] shows the likelihood distribution in e+jets 
and /i+jets data. 

The measurement is consistent with the standard 
model expectation of zero, and is the first such mea- 
surement of the mass difference between a quark and 
its antiquark. 
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Table I Summary of systematic uncertainties on Ami. 



Source 



Uncertainty (GeV) 



Physics modeling 

Signal ±0.85 

PDF uncertainty ±0.26 

Background modeling ±0.03 

Heavy flavor scale factor ±0.07 

6- fragmentation ±0.12 

Detector modeling: 

fa/light response ratio ±0.04 

Jet identification ±0.16 

Jet resolution ±0.39 

Trigger ±0.09 

Overall jet energy scale ±0.08 

Residual jet energy scale ±0.07 

Muon resolution ±0.09 

Wrong charge leptons ±0.07 

Asymmetry in bb response ±0.42 

Method: 

MC calibration ±0.25 

^-tagging efficiency ±0.25 

Multijet contamination ±0.40 

Signal fraction ±0.10 



Total (in quadrature) 



±1.22 



Table II Summary of systematic uncertainties for the lep- 
ton±jets matrix element mass measurement. 



Source 



Uncertainty (GeV) 



Physics modeling: 

Signal modeling ±0.40 

PDF uncertainty ±0.14 

Background modeling ±0.10 

^-fragmentation ±0.03 

Detector Modeling: 

6/light response ratio ±0.83 

Jet identification and resolution ±0.26 

Trigger ±0.19 

Residual jet energy scale ±0.10 

Muon resolution ±0.10 

Method: 

MC calibration ±0.26 

^-tagging efficiency ±0.15 

Multijet contamination ±0.14 

Signal contamination ±0.13 

Signal fraction ±0.09 



Total 



±1.07 



4.1. Performance on Simulated Data, 
Calibration of Method 
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Figure 3: Amt measured in data for a.) e±jets b.) 
mu±jets. 



4. Lepton+Jets Matrix Element Top Mass 
Measurement 



As in the mass difference measurement, ensembles 
of pseudoexperiments are constructed to test the per- 
formance of the mass measurement technique, and bi- 
ases on mt, 5mt, JES, and 8 JES are removed (see 

Fig. m. 
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Figure 4: < mf casurcd > and pull RMS vs mf ue . 



The ME mass measurement uses the same method 
as the mass difference measurement, but with m t — 
mi, and with the reconstructed jet energies adjusted 
by a constant factor JES, which is treated as a vari- 
able in the likelihood, and which calibrates the jet 
energy scale using the hadronic W mass. 



4.2. Systematic Uncertainties 

The largest systematic uncertainty come from un- 
certainties in the calibration of the energy of b-jets. 
Table ITT1 lists the systematic uncertainties. 
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4.3. Measurement in Data 

The measurement in data is @: 

1.0 fb" 1 m t = 171.5 ± lA(stat) ± 1.8(syst) GeV 
2.6 ftT 1 m t = 174.8 ± l.O(stat) ± 1.6(syst) GeV 
3.6 fb" 1 m t = 173.8 ± 0.8(stat) ± 1.6(syst) GeV 

(6) 



The background from fake leptons is controlled with 
lepton quality cuts, the Z/j* — > ee background is con- 
trolled by requiring large Et and removing events 
within the Z peak, the Z/j* — > /i/i background is 
controlled by requiring a large missing ET significance 
(E T /(2a 2 )) where the error is calculated per event us- 
ing the known resolutions of the objects in the event. 
The Z/j* — > tt background is reduced in the efj, chan- 
nel with a cut on the sum of the leading lepton px plus 
the sum of the px of the two leading jets. 



Here the error coming from the in-situ JES has been 
added to the systematic uncertainty. Figure [5] shows 
contours of the likelihood in the M top vs. JES plane. 



6. Dilepton Matrix Element Mass 
Measurement 
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Figure 5: Contour plot showing mt vs JES in data 



The Matrix Element is employed here in the same 
way as for the lepton+jets measurement, the differ- 
ence is that there is no JES factor applied to jets, as 
there is no hadronic W in dilepton events. 

6.1. Performance on Simulated Data, 
Calibration of Method 

As in the 1+jets mass measurement, ensembles of 
pseudoexperiments are constructed to test the perfor- 
mance of the mass measurement technique, and biases 
on mt and Smt are removed (see Fig. 



p1 0.9157±0.05B7 



D0 Run lib prel., L=2.5fb' 



X J ' ™J1 3.909 / 1 

pO 1.142±0.0141 



true m -172.5 (GeV) 



true m, -172.5 (GeV) 



5. Dilepton Event Selection 



Figure 6: < m ™ asurcd > and pull RMS vs m^ rue . 



There are five dilepton channels, three of which 
use only well identified electrons and muons, the ee, 
e[i, and fi^L channels, and two which require one well 
identified electron or muon matched with an isolated 
charged track, the e+track and /i+track channels. 

All five channels require lepton pr > 15 GeV, all 
channels require two jets with px > 20 GeV. 

The main backgrounds are: 

• events with fake leptons (QCD, VF+jets) 

• Z/j* — > ee/n/j, 

• Z/j* — > TT — > l~VlV T l + Vl'V T 

• WW -» l~vil + 'vi' 



6.2. Systematic Uncertainties 

The largest systematic uncertainties come from the 
jet energy scale calibration (both the light jet scale, 
and the b-jet light -jet ratio) (see Fig. [7]). 

6.3. Measurement in Data 

The measurement in data is: 

1.0 fb" 1 m t = 171.7 ± 6A(stat) GeV 
2.6 fb" 1 m t = 176.1 ± 3.9(stat) GeV 
3.6 fb" 1 m t = 174.8 ± 3.3(stat) ± 2.6(syst) GeV 

Figure [8] shows the likelihood curves for the 1 fb" 1 
and 2.6 fb" 1 . 
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Figure 7: Systematic uncertainties for dilepton matrix el- 
ement mass measurement. 



Where fix) and f{x) are parton distribution func- 
tions for the incoming partons, and p(E^\mt) and 
p(E~\mt) give the probability for the top to decay 
into a lepton with energy E\ in the rest frame of the 
top, given by the expression: 



p(E* e \m t ) 



AE^m t (m 2 — - 2E^m t ) 



(mf - ml) 2 - m 2 w {m 2 + mg) 



(9) 

Both methods sum the weights produced at each 
hypothetical top mass for some number of resolution 
smearings ( 100), and for the different solutions and 
jet-lepton assignments (see Fig. [9]): 
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Figure 8: Likelihood distribution in data for Runlla, and 
Runllb. 



7. Dilepton Template-based Mass 
Measurements 

The two template based methods start by using the 
kinematic information in the event, along with a hy- 
pothesized top mass, to reconstruct one or more kine- 
matic tt configurations consistent with the observed 
event. In addition, both methods smear the lepton 
and jet momenta within their known detector resolu- 
tions. 

The neutrino weighting (vWT) method samples the 
neutrino rj distribution from expected distributions, 
then uses this information along with the lepton and 
jet momenta to solve the event. Finally, a weight is 
constructed which compares the fix calculated using 
the neutrino momenta to the measured $t '■ 
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Figure 9: Weight distribution in i/WT. 

The i/WT method takes the mean and RMS of this 
sum-of-weights distribution, and constructs 2D dis- 
tributions of signal and background known as tem- 
plates (see Fig. [TUl note that vWT uses two correlated 
techniques referred to as i/WT^and z/WT/), while the 
MWT method takes the mass at which the weight 
curve is maximized m™ axw , and forms ID templates 
(see Fig. [TTj) . The templates are then used in a likeli- 
hood to measure the top mass. 
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(7) 

The matrix weighting method uses the measured I£t 
and uses an analytical technique to solve directly for 
the kinematics[6|[7[. Each solution is then weighted 
by the formula: 



w = f{x)f{x)p{E* t \m t )p{Ej\m t ) 



(8) 



7.1. Performance on Simulated Data, 
Calibration of Method 

As before, ensembles of pseudoexperiments are con- 
structed to test the performance of the mass mea- 
surement technique, and biases on m t and Suit are 
removed. Table UTT1 gives the slope and offset for these 
calibrations, while Fig. [T2l show the uncalibrated plots 
for the i^WT and MWT. 
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Table III Slope (a) and offset (/?) from the linear fit to the pseudoexperiment results for the 21, ltrk, and combined 
dilepton channel sets using the MWT and i/WT methods. 



Method Channel Slope: a Offset: (3 [GeV] Pull width 



Expected statistical 
uncertainty [GeV] 
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Figure 10: Template distribution in i/WT. 



Figure 12: Calibration curves i/WT and MWT. 



>, 0.2 

"</) 

|0.15 



2 0.1 



CO 

n 
o 



0.05- 



D0, MWT 
m t =170 GeV 



100 120 140 160 180 200 220 240 

Top Mass [GeV] 



7.3. Measurement in Data 

The measurement in data is 

i/WT m t = 176.2 ± 4.8 (stat) ± 2.1(syst) GeV 
MWT m t = 173.2 ± 4.9{stat) ± 2.0(syst) GeV 
Combined m t = 174.7 ± AA(stat) ± 2.0{syst) GeV 

(11) 

Figure [TJ shows -logC for the vWT and MWT 
methods. 



Figure 11: Template distribution in MWT. 



7.2. Systematic Uncertainties 

As in the dilepton matrix element analysis, the 
dominant systematic uncertainties are the jet energy 
calibration and b-to- light jet ratio. Table IIVI lists the 
systematic uncertainties. 



8. Combination of D0 Top Mass 
Measurements, Conclusions 

The D0 experiment has performed the first mea- 
surement of the mass difference between the top quark 
and its antiquark, which is the first such measurement 
for any quark. In addition the D0 experiment has 
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Table IV Summary of systematic uncertainties for the 
combined analysis of all dilepton channels. The i/WTt, 
vWTf, and MWT method results are shown. 



summarizes the current mass measurements for D0. 



Source of uncertainty 


vWT h 


i/WT/ 


MWT 




[GeV] 


[GeV] 


[GeV] 


b fragmentation 


0.4 


0.5 


0.4 


Underlying event modeling 


0.3 


0.1 


0.5 


Extra jets modeling 


0.1 


0.1 


0.3 


Event generator 


0.6 


0.8 


0.5 


PDF variation 


0.2 


0.3 


0.5 


Background template shape 


0.4 


0.3 


0.3 


Jet energy scale (JES) 
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1.6 


1.2 


6/light response ratio 


0.3 


0.4 


0.6 


Sample dependent JES 


0.4 


0.4 


0.1 


Jet energy resolution 


0.1 


0.1 


0.2 


Muon/track pr resolution 


0.1 


0.1 


0.2 


Electron energy resolution 


0.1 


0.2 


0.2 


Jet identification 


0.4 


0.5 


0.5 


MC corrections 


0.2 


0.3 


0.2 


Background yield 


0.0 


0.1 


0.1 


Template statistics 


0.8 


1.0 


0.8 


MC calibration 


0.1 


0.1 


0.1 


Total systematic uncertainty 


2.1 


2.3 


2.0 
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